JOURNAL OF GEOPHYSICAL RESEARCH, VOL. ???, XXXX, DOI: 10.1029/, 


A topside equatorial ionospheric density and 
composition climatology during and after extreme 
solar minimum 

J. H. Klenzing, 1 F. Simoes, 1 S. Ivanov, 2 R. A. Heelis, 3 D. Bilitza, 4 ’ 5 R. F. 

Pfaff, 1 and D. E. Rowland 1 


J. H. Klenzing, Space Weather Lab / Code 674, Goddard Space Flight Center, Greenbelt, 
Maryland, USA (jeffrey.klenzing@nasa.gov) 

1 Space Weather Lab / Code 674, 

Goddard Space Flight Center, Greenbelt, 

Maryland, USA 

department of Physics, Georgia Institute 
of Technology, Atlanta, Georgia, USA 

3 Center for Space Sciences, University of 
Texas at Dallas, Richardson, Texas, USA 

4 Heliophysics Laboratory / Code 672, 

Goddard Space Flight Center, Greenbelt, 

Maryland, USA 

5 Department of Physics, George Mason 
University, Fairfax, Virginia, USA 


DRAFT 


September 6, 2011, 3:01pm 


DRAFT 



X - 2 


KLENZING ET AL.: TOPSIDE EQUATORIAL IONOSPHERIC DENSITY 


4 Abstract. During the recent solar minimum, solar activity reached the 

5 lowest levels observed during the space age. This extremely low solar activ- 

e ity has accompanied a number of unexpected observations in the Earth’s iono- 
7 sphere and thermosphere when compared to previous solar minima. Among 
s these are the fact that the ionosphere is significantly contracted beyond ex- 

9 pectations based on empirical models. Climatological altitude profiles of ion 

10 density and composition measurements near the magnetic dip equator are 

11 constructed from the C/NOFS satellite to characterize the shape of the top- 

12 side ionosphere during the recent solar minimum and into the new solar cy- 

13 cle. The variation of the profiles with respect to local time, season, and so- 

14 lar activity are compared to the IRI-2007 model. Building on initial results 

15 reported by Heelis et al. [2009] , here we describe the extent of the contracted 

16 ionosphere, which is found to persist throughout 2009. The shape of the iono- 

17 sphere during 2010 is found to be consistent with observations from previ- 

18 ous solar minima. 
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1. Introduction 

19 The solar minimum between cycles 23 and 24 has been an unusual period of solar activ- 

20 ity. The minimum was expected to occur in 2006, but instead solar activity continually 

21 decreased throughout 2007 and 2008 [Russell et a/., 2010]. Traditional proxies for solar ac- 

22 tivity such F10.7 (the flux of solar radiation at 10.7 cm wavelength) have “bottomed out,” 

23 while actual measurements of the EUV flux have continued to decrease [Araujo- Pradere 

24 et al , 2011; Chen et a/., 2011]. Some long-term climate modelers have even speculated 

25 that the deepest part of this minimum could be used to better understand the Maunder 

26 minimum [Schrijver et a/., 2011]. 

2? A number of surprising observations in the ionosphere and thermosphere has accom- 

28 panied this period of extremely low solar activity. The thermospheric density was found 

29 to reach record lows based on the analysis of the orbital decay of numerous satellites 
so [ Emmert et al. , 2010] and by in situ measurement of the neutral scale-height [Haaser et 
si a/., 2010]. Solomon et al. [2011] showed that this reduction in thermospheric density was 

32 largely due to low solar activity and that other secular variations (such as geomagnetic 

33 activity) were small in comparison. 

34 Heelis et al. [2009] showed that the ionosphere was contracted as well, with the transition 

35 height between H + and 0 + being significantly lower than predicted by the IRI model, 
se Additionally, the topside nighttime ion temperatures have been found to be relatively 
s? cold compared to IRI at an altitude of 400 km, as low as 600 K [Coley et a/., 2010]. The 
ss average ion drift in the topside ionosphere has been found to be significantly different 
39 from previous observations, including other solar minima. The ExB drift climatology 
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40 observed by the C/NOFS satellite was found to differ from the Fejer-Scherliess model 

4 1 [e.g., Scherliess and Fejer , 1999], including downward afternoon drifts in some regions as 

42 well as a weak to non-existent pre-reversal enhancement [Pfaff et al , 2010]. 

43 The behavior of the topside ionosphere during solar minimum has been well- documented 

44 through in situ measurements [Greenspan et al. , 1994; West et al. , 1997], topside sounders 

45 [Benson and Bilitza , 2009], and ground-based radar [Hysell et al. , 2009]. However, because 

46 solar activity is lower during the cycle 23/24 minimum than during the last few solar 

47 cycles, current empirical models must extrapolate based on previous observations. Long- 

48 term monitoring of ionospheric density by the CHAMP and GRACE satellites reveal that 

49 the IRI-2007 model overestimates the expected density leading up to and including the 
so recent solar minimum [Liihr and Xiong , 2010]. 

5 1 In this study, the ion density and composition data from the C/NOFS satellite are used 

52 to construct climatological maps of density and composition as a function of altitude and 

53 local time near the magnetic dip equator. The maps are divided into season and solar 

54 activity in order to understand the shape of the topside ionosphere during this extreme 

55 solar minimum, as well as its evolution on the journey back to solar maximum. These 
se climatology maps are compared to the results from the IRI-2007 model. The highly 

57 contracted ionosphere is found to persist throughout 2009, well into the new solar cycle. 

2. Measurements and Models 

58 The Communication/Navigation Outage Forecast System (C/NOFS) satellite is part of 

59 a space weather mission led by the US Air Force Research Laboratory to locate, under- 
eo stand, and predict equatorial ionospheric scintillation [de La Beaujardiere et al ., 2004]. 
6 i The C/NOFS satellite was launched in April 2008 into a 13° inclination orbit with perigee 
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62 near 400 km and apogee near 860 km. This elliptical orbit allows for a sampling of ion 
es density over multiple scale heights of the topside ionosphere. The C/NOFS perigee pre- 
64 cesses through all solar local times roughly once every 65 days. C/NOFS is equipped with 
es multiple instrument suites designed to study the ion and neutral populations and their 
ee effect on the propagation of communication signals. 

e? This study will focus on the total ion density (electron density for a quasi-neutral 
es plasma) and the H + and 0 + components. The total density and composition are obtained 
eg from a Retarding Potential Analyzer (RPA), part of the Coupled Ion-Neutral Dynamics 

70 Investigation (CINDI) suite of instruments on board C/NOFS. The well-established RPA 

71 technique consists of using a series of biased grids to select certain energies of ions to 

72 measure as a current [Heelis and Hanson , 1998]. By sweeping over a range of voltages, 

73 the relative contribution of each ion species (along with ion drift velocity and temperature) 

74 can be calculated. 

75 The International Reference Ionosphere (IRI) is considered the international standard 

76 model for calculating empirically- derived ionospheric parameters based on both ground- 

77 based and satellite measurements [Bilitza and Reinisch , 2008]. IRI was founded as a joint 

78 project between the Committee on Space Research (COSPAR) and by the International 

79 Union of Radio Science (URSI). An empirical model was chosen for comparison to the 
so C/NOFS observations in order to better illustrate the differences between the topside 
si density variations when compared to previous solar minima. The IRI model can generate 

82 estimated values of density and ion composition for a given input of solar activity, which 

83 is described by the geophysical indices Rz (based on the sunspot number) and IG (based 

84 on the ionospheric response) [ Bilitza , 2000] . 
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ss The international sunspot number (Rz, also referred to as Ri after the main observation 
se platform moved from Zurich to Brussels in 1980) is a weighted average of the number of 

87 sunspots observed on the surface of the sun [Clette et a/., 2007]. The weighting is designed 

88 to account for the differences between individual sunspots and clusters of sunspots. The 
sq IG index was developed by Liu et al. [1983] to provide an estimate of the peak F 2 -region 

90 density based on the International Radio Consultative Committee (CCIR) maps of the 

9 1 ionosphere. To do this, the index uses a weighted average of ionosonde measurements 

92 around the world similar to the method introduced by Minnis and Bazzard [I960]. This 

93 index is scaled to produce a “global effective sunspot number.” However, since this cal- 

94 culation is based on the actual measurements of the ionosphere, it does not have a lower 

95 limit (unlike Rz, which by definition cannot be less than zero). During the prolonged solar 

96 minimum between cycles 23 and 24, IG is often negative. The IG index is provided by 

97 the UK World Data Center and is available in either monthly averages (IG) or 12-month 

98 averages (IG 12 ). 

99 The NeQuick topside model is used to generate expected values of ion density and 

100 composition every five seconds along the orbit track of C/NOFS for every orbit through 

101 the end of 2010. The NeQuick model was chosen due to its excellent performance when 

102 compared to topside sounder measurements from the Alouette and ISIS missions relative 

103 to the other topside options included in IRI [Bilitza et a/., 2006; Bilitza , 2009]. For this 

104 study, the 12 -month running averages of Rz (referred to as Rz i2 ) and IG (IG 12 ) are used. 

105 Only dates with the definitive values of Rz 12 and IGi 2 (through January 2011 at the time 

106 of this writing) were used. 


DRAFT 


September 6, 2011, 3:01pm 


DRAFT 



KLENZING ET AL.: TOPSIDE EQUATORIAL IONOSPHERIC DENSITY 


X- 7 


107 To illustrate the prolonged nature of the recent solar minimum, these activity proxies 
ice are compared for the last three solar cycles in Figure 1. Panel (a) shows the values of 
loo F10.7A (the 81-day average of F10.7) for 36 months around the minima between cycles 
no 20 and 21 (1975, plotted in blue), cycles 21 and 22 (1986, plotted in green), cycles 22 and 
in 23 (1996, plotted in orange), and cycles 23 and 24 (2008, plotted in red). The values of 
n 2 Rzi 2 and IG 12 for the same three periods are plotted in panels (b) and (c), respectively, 
ns All three indices show that solar activity is deeper and longer than in previous years. 

114 A sample portion of the data used in this study (~3 consecutive orbits) is shown in 
ns Figure 2. The CINDI measurements are shown as solid lines, the predicted IRI-2007 values 
no are shown as dashed lines. Two examples are given to illustrate the effect of precession 
117 on the density variations: the 17 Nov event (a) shows three consecutive orbits during a 
ns period when perigee is near local noon; the 19 Dec event (b) shows a similar section of 

119 data when perigee is near local midnight. Note that the CINDI measurements are more 

120 variable than the IRI predictions. 

3. Technique for Reconstructing Topside Profiles 

121 The density measurements from C/NOFS are averaged together to create climatological 

122 altitude profiles. Because the elliptical orbit precesses through a variety of longitudes, 

123 the reconstructed profiles cannot be thought of as the ionospheric profile at any given 

124 location. Rather, these are average profiles that neglect longitudinal variations, tidal 

125 effects and local magnetic anomalies. (Such effects are small relative to the altitude and 

126 local time effects and will be the topic of a future study.) In order to be certain that any 

127 comparisons between the data and model are on equal footing, the algorithms used to 
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128 reconstruct average topside profiles from the CINDI data are also used on the expected 

129 density values generated by IRI along each orbit of C/NOFS as shown in Figure 2. 

iso The C/NOFS satellite undergoes a complete precession of perigee through all local 

131 times roughly once every 65 days. In order to better smooth out the variations due to 

132 the daily longitudinal precession, a period of 91 days is used for the reconstructed topside 

133 data. Data that is noisy or contains localized features such as plasma density depletions or 

134 enhancements are removed from the averages in order to approximate a true background 

135 density. To accomplish this, a Savitsky-Golay filter [Savitsky and Golay , 1964] is used on 
ise the C/NOFS data to determine the smoothness of the dataset. The smoothing window is 
is? 241 points wide (containing roughly two minutes of data), and the smoothing function is a 
ise third-degree polynomial. Because the filtered values represent an average of the perturbed 
i39 and background densities, filtered values that differ by more than 0.5% from the measured 
uo value are removed from the averages. 

i4i The profiles are calculated every 0.5 hours in local time with a 10 km resolution. To 
U2 smooth out the data, there is a significant overlap between each bin. (The bins are 2.5 
us hours wide and 50 km high.) Only points within ±2.5 deg magnetic dip latitude are used 
144 in the reconstruction, but all longitudes are used in the averages. To remove the effects 
us of magnetic substorms, only quiet times where Kp < 3 are used. 

U6 Figure 3a shows a sample dayside profile for the December Solstice of 2008 as recon- 
147 structed from CINDI data (shown in black) and from the IRI predictions (shown in green) . 
us The dashed lines represent the first and third quartiles for each bin. Part of this vari- 

149 aiton in the data set is due to the longitudinal variations. A sample post-sunset profile is 

150 shown in Figure 3b. Note that there is a sharp in the vertical gradient for the measured 
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isi profile, while the IRI prediction varies smoothly. A similar effect was observed for higher 

152 magnetic latitudes in Sibanda and McKinnell [2011]. 

153 Figure 4 shows the average composition associated with the profiles in Figure 3. The 

154 solid lines represent the measured components of H + (red) and 0 + (blue), and the dashed 

155 lines represent the IRI expectations. Note that for both the dayside and night side profiles, 

156 the measured concentration of 0 + is consistently lower than the expected value. This is 

157 consistent with previous findings that the ionosphere is contracted more than expected 

158 in the recent solar minima [Heelis et al, 2009; Liihr and Xiong , 2010]. However, the H + 

159 component may be either larger or smaller than expected; it is consistently larger on the 
leo dayside profile in the altitude range of the C/NOFS satellite. 

lei The transition height between 0 + and H + can be inferred from the figure by noting 
i 62 where the red curve crosses the blue curve for a given profile. The transition height for 
les the CINDI profile in Figure 4b is is ~50 km lower than the expected value based on IRI. 

164 This is consistent with the findings by Heelis et al [2009] that the transition height is 

165 lower than expected during the recent solar minimum. 

lee Figure 5 shows the composition profiles for the December solstice of 2010, two years 
16 ? after the deepest part of the recent solar minimum. The estimates of both components 
les match the IRI predictions much better for the dayside profiles (5a), as well as the 0 + 
i69 component for the post-sunset profile (5b). The H + component is low by a factor of 4 for 
no the upper altitudes. 

171 Density and composition profiles as shown in Figures 3-4 are generated for every 0.5 

172 hours of local time for 91-day seasons ranging from the December solstice of 2008 to the 

173 December solstice of 2010. The variations of the topside density profiles with respect to 
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174 local time and season are discussed in the following two sections. The seasonal divisions 

175 and the associated average solar indices are listed in Table 1. (The final season is cut 

176 short by 4 days due to the current availability of definitive values for the Rz i2 and IGi 2 

177 indices for driving the IRI model). The two equinoctial seasons remain separate in order 
ns to better capture the effects of the slowly increasing solar activity. 

4. Variation with Local Time 

179 Figure 6 is the summary plot for the December Solstice of 2008. (This corresponds 
iso to the deepest part of the cycle 23/24 minimum.) Panel (a) shows the average total 

181 density as measured by CINDI as a function of altitude and solar local time, and panel (b) 

182 shows the equivalent average from the IRI-2007 values generated over the C/NOFS orbits. 

183 The transition height between H + and 0 + as calculated from the average reconstructed 

184 profiles is plotted over the contour maps as a solid (CINDI) or a dashed (IRI) black line. 

185 Figure 6c is the ratio of the IRI average density to that computed from CINDI. (Both 

186 the measured and modeled transition heights are included in this panel for reference.) 
is? Similarly reconstructed profiles for the concentrations of 0 + and H + are shown in Figures 

188 6d-i. 

189 Note that over the full range of altitudes and local times covered by the C/NOFS 

190 satellite, there can be found regions where the IRI-2007 model will either overestimate 

191 or underestimate density, whether it be total density (panel c), 0 + concentration (f), 

192 or H + concentration (i). Liihr and Xiong [2010] found that IRI tended to overestimate 

193 total density during solar min; this study was conducted with the CHAMP and GRACE 

194 satellites, which are in circular orbits at 310 and 490 km altitude, respectively. Similarly, 

195 for a fixed local altitude near the C/NOFS perigee, Figure 6c shows that the IRI model 
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196 overestimates total density for all local times except near the dawnside terminator. Av- 

197 eraging over all local times at 490 km, IRI overestimates the C/NOFS density by about 

198 80% for this time period, which is consistent with the GRACE results. 

199 Several additional features of the contracted ionosphere are clearly seen in Figure 6. In 

200 particular, the post-sunset electron density is lower than predicted by up to a factor of 4 

201 (6c), and the concentration of 0 + is generally smaller than predicted for all local times 

202 (6f), except near the dawnside terminator. The measured concentration of H + is larger 

203 than estimated by IRI on the dayside profiles and for the nightside below the transition 

204 height. 

205 Figures 7 and 8 show similar density maps for the December solstices of 2009 and 2010, 

206 respectively. It is readily apparent that the densities for the March equinox of 2010 are 

207 much closer to the expected values. For instance, both the ratio plots for total density (c) 

208 and 0 + concentration (e) are significantly closer to one when compared to the previous 

209 year. However, there are still some discrepancies between model and data, such as near 

210 sunrise and above the transition height post sunset. The estimates of H + for 2010 are 

211 high on the nightside and low on the dayside. 

5. Variation with Season and Solar Activity 

212 Climatological maps similar to those shown in Figures 6-8 were generated for each 

213 season shown in Table 1. In order to better illustrate the effects of seasons and solar 

214 activity, certain metrics will be adopted. Figure 9 shows the vertical “total electron 

215 content” (TEC) between 400 and 800 km. (This should not be confused with the total 

216 electron content in the typical sense, since we only observe over a relatively small range of 

217 altitudes above the F2-peak. However, it is instructive to display the relative changes in 
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218 this metric.) Because the height of apogee drops over the lifetime of the satellite (822 km 

219 as of 31 January 2011), 800 km was chosen as the upper limit for the integrated density 

220 measurements for the purposes of this metric (although the total density above 800 km is 

221 minimal). The integrated vertical density is plotted in TEC units (TECU), or units of 10 16 

222 m -3 . Each panel represents a given season and contains the TEC400-800 calculated from 

223 the CINDI measurements for each available year plotted as solid lines, with data from 

224 2008 plotted in green, 2009 in purple, and 2010 in orange. The corresponding TEC400-800 

225 as calculated from IRI-2007 is shown in dashed lines. 

226 The measured densities are significantly lower for all local times through the December 

227 solstice of 2009. For the March equinox of 2010 (the orange line in Figure 9c), the 

228 integrated density is ~25% larger than the model on the dayside. For the two following 

229 seasons in 2010, the modeled values better approximate the CINDI measurements than in 

230 the previous seasons. Additionally, the CINDI measurements in the December solstice of 

231 2009 are still very close to those from 2008, while there is a dramatic increase in density for 

232 the March equinox of 2010. This corresponds to a rise in F10.7 above 80 sfu. Additionally, 

233 the average effective sunspot number (IG12) based on the ionospheric activity more closely 

234 matches the measured sunspot number (Rz i2 ) for this period. 

235 Another interesting feature is that the measured ion density in March equinox of 2010 

236 is larger than the corresponding density in the September equinox of the same year. The 

237 TEC plots from these two periods have been replotted in a single panel in Figure 10 to 

238 better illustrate this asymmetry. All three of the solar /ionospheric activity proxies are 

239 larger in for the September equinox than the March period (see Table 1). Accordingly, 

240 the IRI-2007 modeled values predict that the densities in September would be larger (the 
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24 1 dashed orange lines) than those from March (dashed purple lines) . This is clearly not the 

242 case in the measured densities (represented by the solid lines) . This equinoctial asymmetry 

243 is similar to that noted in the COSMIC TEC data by Liu et al. [ 2010 ]. A similar effect 

244 was recently reported in the vertical drift data from the ROCSAT -1 satellite [Ren et al . , 

245 2011 ]. However, we should remember that the ROCSAT data is from a period when 

246 the solar activity was much higher, and an investigation into this equinoctial asymmetry 

247 utilizing the C/NOFS vertical drift data will be required to fully understand this ion 

248 density asymmetry. 

249 The transition height between H + and 0 + is calculated from the reconstructed profiles 

250 for both data and model. These are shown with respect to seasonal and temporal vari- 

25 1 ations in Figure 11 (similar to the integrated densities presented previously). Note that 

252 there is no significant difference between the nightside transition height as predicted by 

253 IRI over the course of the mission (the dayside transition heights are typically outside 

254 of the range of the C/NOFS satellite). The dayside transition heights for 2010 are still 

255 lower than predicted by the models, but it should be noted that these are consistent with 

256 observations from Atmospheric Explorer from a previous solar minimum [Gonzalez et a/., 

257 1992 ]. The nightside transition heights for 2010 are consistent with that predicted by IRI. 

6. Discussion 

258 The reduced densities observed during the extreme portion of the recent solar minimum 

259 could be explained by any combination of the following effects: 

260 1 . The height of the F 2 layer (hmF 2 ) is lower than predicted. 

261 2 . The density of the F 2 peak (NmF 2 ) is lower than predicted. 
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262 3. The shape of the topside ionosphere is different. 

263 Because the peak of the F layer is below the C/NOFS perigee (400 km) for most of 

264 the mission, we cannot comment on the cause from this data alone. (C/NOFS travelled 

265 below the F peak for the first time in April 2011.) Figure 12 is provided to illustrate this 

266 problem. An initial ion density profile is generated from IRI-2007 (shown in black) for 

267 the December solstice of 2008 for 1400 local time. Two modified profiles are created that 

268 would lead to the observed ~67% overestimate in TEC400-800 shown in Figure 9. The 

269 first is created by simply scaling the density by a factor of 0.6 (shown in blue), the second 

270 is created by moving the F peak down by 65 km (shown in red). 

271 While recent studies using ionosonde data have shown that NmF2 reached record low 

272 measurements during the recent solar minimum [Liu et al , 2011], IRI-2007 is found to 

273 predict this density very well, with the standard deviations being comparable to previous 

274 solar cycles [Bilitza et al. , 2011]. This is due to the fact that the IRI model predicts the 

275 peak density based on the IG12 index, which is in itself a global average of ionosonde 

276 measurements. 

277 The topside ionospheric density is controlled not just by solar radiation, but by a balance 

278 of chemical and dynamic processes. In the topside ionosphere, the creation of H + ions is 

279 primarily due to charge exchange with 0 + [Rishbeth and Garriott , 1969]. 

0 + + H^O + H + ( 1 ) 

280 To first order approximation, the relation between the the ion components will then 

281 depend on the densities of the neutral components. 
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< 2 > 

282 The factor of 9/8 is due to statistical differences in the forward and reverse reaction 

283 rates. The increased concentration of H + in the topside ionosphere is consistent with the 

284 increased ratio of neutral [A/]/[0] observed in the upper thermosphere [ Haaser et al. , 2010] . 

285 A recent study by Hysell et al. [2009] compared topside profiles from the Jicamarca Radar 

286 Observatory with the SAMI2 model and concluded that the shape of the H + fraction is 
28 ? also affected by the ExB drift time history as well. 

288 This decreased post-sunset density may be partially related to the altered vertical drift 

289 climatology during the recent solar min. Unlike the Fejer-Schierless model, the vertical 

290 drift is found to be downward in the afternoon, and the large upward drift around sunset 

291 known as the pre-reversal enhancement is largely absent in the 2008 and 2009 data [Pfaff 

292 et a/., 2010]. The density structure may also be due to different climatologies in the 

293 meridional winds. 

7. Summary and Conclusions 

294 A statistical study of the ion density and composition in the topside ionosphere near 

295 the magnetic dip equator during the recent solar minimum was conducted. The major 

296 findings are the following: 

297 1. While the overall ionosphere was found to be contracted relative to empirical ex- 

298 pectations, the ratio of the expected density to the measured density was found to be a 

299 strong function of altitude and local time, including some areas (such as ~800 km just 
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3oo before dawn) where the average measured ion density was higher than predicted by as 
3d much as a factor of four. 

302 2 . During this contracted phase, [H + ] is found to be greater than predicted by IRI -2007 

303 for all observed altitudes (400 to 850 km) on the dayside and below the transition height 

304 for the night side. 

305 3 . The shape of the topside nighttime ionosphere between 400 and 850 km was found 

306 to be different from the predicted shape. The profile generated by IRI varies smoothly, 

307 while the data shows a sharp change in the vertical gradient associated with the lower 

308 transition height. 

309 4 . The post-sunset ion density decreased more rapidly than expected based on previous 

310 solar minima. This may be related to a different drift climatology observed with the 

311 C/NOFS satellite during extreme solar min as previously reported by Pfaff et al. [ 2010 ]. 

312 5 . This highly contracted ionosphere persisted until the March equinox of 2010 , over a 

313 year into the new solar cycle. The transition heights observed in 2010 are are consistent 
3u with observations from previous solar minima. 

sis 6. The geophysical indices used to drive the IRI model, Rz i2 and IG12, are both signif- 

316 icantly lower than in previous solar minima. The previously reported tendency of IRI to 

317 overestimate density during the extreme solar min is not a deficiency of the chosen input 
sis indices, but rather illustrates the fact that we have not observed the ionosphere during 

319 such a low period of solar activity. The reconstructed topside profiles from C/NOFS can 

320 be used as an additional constraint on future versions of IRI. 

321 The C/NOFS satellite provides a unique look at the shape of the topside ionosphere. 

322 The topside data from C/NOFS during this unprecedented low in solar activity could 
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323 be used as a constraint on future empirical models. Future studies will include regional 

324 case studies for comparison with ground-based measurements, as well as variations with 

325 longitude and magnetic latitude. Additionally, the reconstructed profiles can be used to 

326 discuss transport phenomena in the topside ionosphere in conjunction with drift clima- 

327 tologies using physics-based models to quantify the relative effects of altered transport 

328 and chemistry during extreme solar minimum. 
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Table 1 . Seasonal divisions for the topside profile reconstructions, including the average solar 


activity represented by F10.7, Rzi 2 , and IG 12 for each period. 


Period 

Range of Days 

F10.7 

RZ 12 

IG 12 

Sep Equinox 2008 

8 Aug 2008 

- 6 Nov 2008 

67.8 

2.3 

-10.1 

Dec Solstice 2008 

6 Nov 2008 

- 4 Feb 2009 

67.3 

2.0 

-9.2 

Mar Equinox 2009 

3 Feb 2009 - 

- 4 May 2009 

69.2 

2.0 

-8.4 

Jun Solstice 2009 

7 May 2009 

- 5 Aug 2009 

70.9 

3.2 

-6.3 

Sep Equinox 2009 

8 Aug 2009 

- 6 Nov 2009 

70.9 

6.2 

-0.8 

Dec Solstice 2009 

6 Nov 2009 

- 4 Feb 2010 

75.5 

8.5 

3.9 

Mar Equinox 2010 

3 Feb 2010 - 

- 4 May 2010 

80.3 

12.5 

10.3 

Jun Solstice 2010 

7 May 2010 

- 5 Aug 2010 

77.9 

16.8 

15.7 

Sep Equinox 2010 

8 Aug 2010 

- 6 Nov 2010 

81.5 

21.0 

19.0 

Dec Solstice 2010 

6 Nov 2010 - 

- 31 Jan 2011 

81.6 

29.1 

26.2 
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Figure 1 . Solar activity near the solar minima for the last three cycles, inlcuding (a) F10.7A, 
(b) Rzi 2 , and (c) IGi 2 . The F10.7A values are 81-day averages, and the Rz i2 and IGi 2 indices 
are 12-month averages. Rz i2 and IGi 2 are used to drive the IRI-2007 model in this study. 
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Figure 2. The variability of C/NOFS density data, along with the expected values based on 
IRI-2007. The top panel shows five hours (roughly three orbits) from 17 Nov 2008 (when perigee 
is at local noon), and the bottom panel shows the same for 19 Dec 2008 (when perigee is at local 
midnight) 
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Daytime (SLT = 14.25) 



Nighttime (SLT = 21 .25) 



Figure 3. A sample reconstructed altitude profile based on the average C/NOFS CINDI data 
(black) for the topside equatorial ionosphere for the December solstice of 2008, along with the 
associated IRI-2007 profile (green). The solid lines represent the median density profile; the 
dashed lines represent the first and third quart iles. The left panel is a dayside profile, centered 
around 14.25 local time; and the right panel is a nightside profile, centered around 21.25 local 
time. 
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Daytime (SLT = 14.25) 



Nighttime (SLT = 21 .25) 



Figure 4. The composition profiles associated with the density profiles from Figure 3. The two 
major component ions are H + (red) and 0 + (blue). Note that while the total ion density matched 
IRI quite well for the dayside profiles, the composition is quite different. For both dayside and 
nightside, the transition height between H + and 0 + is lower than predicted. 
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Daytime (SLT = 14.25) 


Nighttime (SLT = 21 .25) 



Figure 5. The same as Figure 4, but for the December solstice of 2010 (2 years later). 
Note that the nightside transition height between 0 + and H + is now very similar for both the 
measurements and the model. 
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(a) September Equinox (b) December Solstice 



Figure 9. The “total” electron content between 400 and 800 km as a function of solar local time. 
These plots capture the seasonal and temporal variation of density for (a) September equinox, (b) 
December solstice, (c) March equinox, and (d) June solstice. Note that the measured densities 
are significantly lower than predicted by IRI until the March equinox of 2010. 
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Figure 10. Selected data from Figure 9 replotted to illustrate the equinoctial asymmetry 
during 2010. The March equinox TEC is larger than the corresponding data in the September 
Equinox after ~10.5 SLT. 
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(a) September Equinox (b) December Solstice 




(c) March Equinox 


(d) June Solstice 



SLT (hr) SLT (hr) 


Figure 11. The transition height between H + and 0 + as a function of solar local time. These 


plots capture the seasonal and temporal variation for (a) September equinox, (b) December 
solstice, (c) March equinox, and (d) June solstice. Note that the daytime transition height as 
predicted by IRI is above the range of the C/NOFS satellite. 
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Figure 12. The effects of changing the position of the F — peak on the observed topside profile. 
An initial profile (black) is generated using IRI-2007. Two altered profiles are included: one 
where NmF2 is scaled down by 60% (blue), and one where hmF2 is moved down by 65 km. The 
apogee and perigee of the C/NOFS satellite are shown as dashed lines. 
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